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Abstract 

To study the convergence of the loop expansion at the high-temperature electroweak 
phase transition, we calculate the 2-loop effective potential of the 3d SU(2)-Higgs 
model in a general covariant gauge. We find that the loop expansion definitely breaks 
down for large £, but converges rather well for smaller values, deep in the broken phase. 
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Recently, both perturbative |], |], || and non-perturbative ||, || studies of the 
high-T electroweak phase transition have been performed, with the indication that 
perturbation theory works poorly in the symmetric phase, due to non-perturbative 
effects. The purpose of this note is to investigate the convergence of loop expansion 
in the broken phase. The method is to calculate the high-T asymptotic of the 2-loop 
effective potential in a general covariant gauge, and to extract from it a quantity which 
is gauge-independent in the full theory. The degree of gauge-dependence at the 2- 
loop level is then expected to tell something about the convergence of the expansion. 
In accordance with this philosophy, we do not convert the loop expansion into an 
expansion in the coupling constants at any stage. The convergence of loop expansion 
in the broken phase has also been studied in ref. ||, by comparing lattice data and 
perturbation theory. 

Our starting point is the action 



S = Jd 3 x \f*F* + (A$) t (A$) + [™ 2 (/i) + 6m 2 ]&<S> + A($ t $) 



1) 



with F§ = diA^-djAf + ge^A^ and A$ = (d i -igr a A'}/2)^. Here the r a :s are the 
Pauli matrices. The only ^-dependent infinity in this theory is removed by 5m 2 , and 
the coupling constants are RG-invariant. Gauge fixing and compensation is achieved 
by adding to eq. [I| the term 



S g = I d 6 x 



1 (<9,4 a ) 2 + diC a d iC a + ge^d^ A\c c 



2i 



(2) 



All the fields have dimension [GeV] 1 / 2 , and A and g 2 have dimension [GeV]. 

The action S + S g is an effective field theory in the sense of ref. , obtained from 
the gauge-Higgs sector of the standard model after integrating out the heaviest degrees 
of freedom, namely, the non-zero Matsubara frequencies and the Aq fieldQ 0, |3|. The 
relations between the parameters of our 3d theory and those of the Standard Model are 
given in ref. |§. To the assumed accuracy, these relations are gauge-invariant H, 0. 
The temperature dependence of coupling constants is A, g 2 oc T, and that of m 2 (^) is 



m\ii) = 7 (T 2 - T 2 ) + -4/ 2m log ( — ) . (3) 

Here /2m is the coefficient of the 2-loop mass-counterterm, proportional to T 2 when 
expressed in terms of 4d couplings. Notice that the gauge fixing parameter £ gets 
renormalized in the reduction step. 

To calculate the effective potential V(<p), one writes $ = [0 3 +i0 4 , y9+0 1 +i0 2 ] r / \/2 
in the action S + S g and neglects terms linear in quantum fields M . This defines a new 



1 We calculated V(<p) also in the theory from which Ag has not been integrated out, but for 
simplicity we discuss only eq. |l| in detail. The differences between the results are addressed below. 
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theory with the masses m\ 
The non-trivial propagators are 



m 2 (fj,) + 3A<^ 2 , m\ = m 2 (/i) + Ayj 2 , and m 2 T 



4 tt (-*M5(*) 



< t> G (-k)<t> G (k) 

Mk)AK~k) 
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ab 



5ij kikj/k ^ ^kikj 
k 2 + m\ 
k 2 + (m^ 
£; 2 (&; 2 + m 2 ) + £m^m 2 ' 

-^k)A 2 {-k) = -<j> 4 (k)Al(-k) 



k 2 + ml 



k 2 k 2 (k 2 + m?,) + ^m\m\ 



G 



2,3,4 



(4) 



i^m^ki 



k 2 {k 2 + m 2 ) + ^m\m\ 



It is seen that for positive £, m 2 must be positive. Hence the calculation is valid only 
for larger values of if than the tree-level minimum, where m\ = 0. 

At this point, our choice of gauge deserves some explanation. In the literature, one 
sometimes uses an i?^-type of a gauge, since then the non-diagonal terms between the 
longitudinal parts of the A^-fields and the Goldstone bosons are cancelled, yielding 
much simpler propagators. This means that one chooses a different gauge for each 
different value of (p in the effective potential. It is not a priori clear that this is allowed 
(for more criticism, see ref. However, at least in the case of scalar electrodynamics, 
it can be proved that the error induced is of higher order in h [TI| (for a more pessimistic 
conclusion, see ref. fill). To be on the safe side, we chose the (l/2£)(<9jA") 2 -gauge for 
our non-Abelian calculation. 

The 2-loop graphs to be calculated are shown in fig. 1. Due to the fact that the 
000- and <f)A A- vertices include at least one Higgs field 0i, and only the Goldstone 
fields can transform into longitudinal vector bosons, four possible graphs (which are 
not shown) vanish. Also, the graph (f2.a) vanishes due to antisymmetry in momentum 
integration. The method of calculation for the remaining graphs is to write 



k 2 {k 2 + m 2 ,) + ^m\m\ 



m\{K. 



+ 



R 2 



1 



k 2 + m\R 2 _ k 2 + m\R\ 



(5) 



where R 2 ± = l/2±yl/4 — £(m T /m 2 ) 2 , and then to use standard Landau-gauge values 
of integrals, as presented e.g. in ref. @, App.B.2] (we use dimensional regularisation). 
Note that R\ can be complex, but this does not matter, since the poles are off the 
integration path and the final result is real. The 1-loop effective potential is 



VUp) = -m 2 (aW + -\(p A - — 
2 vwy 4 y 12tt 



[6m 3 T + m? + 3ml(R 3 + + R 3 _)} 



(6) 



Due to lack of space, the 2-loop part is presented explicitely elsewhere (ref. |]1J 

There are two basic ways to check the 2-loop result. First, in the limit £ — > 
it reproduces the Landau-gauge result as given in ref. ||. Second, due to the facts 
that we are using a gauge-invariant regularization procedure, and the term $t$ is 
gauge-invariant, the counterterm 5m 2 must be gauge-independent. Indeed, the terms 
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proportional to g 4 £ and g 4 £ 2 from graphs (a), (fl), (dl), (x3), (x4) and (x5) cancel, 
and the divergent piece is 



H~ 4e |> 2 [51 



16tt 2 4e 



<7 4 + 9A<7 2 -12A s 
16 



+ 3g 2 m 2 (fi) \ . (7) 



The quantity inside the square brackets is just / 2m appearing in eq. |3|. 

As was noted above, the integrals appearing in the calculation of V(ip) are defined 
only for m\ > 0, if £ 7^ 0. Since we are interested in making calculations even below 
the temperature where the symmetric phase is metastable, and in our 3d theory this 
roughly means that m 2 (/i) < 0, there is a region of (p where raf, = m 2 (/i) + Xip 2 — > 0. It 
is seen from eq. |] that in this limit integrals including Goldstone bosons become IR- 
divergent, with a gauge-dependent coefficient. In the effective potential, this causes 
part of the diagrams to diverge, as m 2 — > 0. The most severe divergences, proportional 
to g 2 (mfym 2 )£ 5//2 , come from diagrams (fl) and (x4), but these cancel. Logarithmic 
divergences proportional to g 2 m\^ 2 log(m 2 ), coming from the same diagrams, also can- 
cel, as do all the divergences proportional to A(mr/m 2 )^ 2 (my/mi)^ 4 , coming from 
diagrams (c), (fl), and (x2). Still, a divergent piece 



lQn 2 



^ 3/2 0-^ 3/4 e) 1/2 ( i+4 ^ 



A^ 
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from diagrams (h4), (f2), and (c) remains. Even the 1-loop potential is non-analytic at 
the m 2 — > limit, since its m 2 -dependence is of the form ^ 3 / 4 m^ 2 my 2 . The divergence 
in eq. [| is gauge-dependent, and would not appear at all e.g. in the i?^-gauge. 



Let us briefly note that in ref. [R|, a gauge-invariant effective potential of the 
composite operator $t$ is calculated at 1-loop level using the saddle point method. 
One of the two saddle points is at m 2 = 0. A 2-loop calculation of W[J] has exactly 
the same Feynman rules and propagators as our calculation, but the set of diagrams 
is different, since in their calculation all the connected vacuum diagrams are to be 
included, whereas only one-particle-irreducible diagrams enter our calculation [0]. If 
their 2-loop calculation is to be successful in the covariant gauge, the divergent piece 
of eq. [| must therefore be cancelled by the extra diagrams. 

For the present calculation, the crucial question is whether the IR- divergences in 
eq. |] completely spoil the perturbation theory even for m| > 0. This cannot be quite 
so since the £ — > 0-limit is the well-behaved Landau gauge. From the formulas for 
I$±, it is seen that one is effectively at this limit when £ <C {m^jlvnT) 2 ~ ^1 9 2 ■ This 
would indicate that only extremely small values of £ are allowed, since typically at the 
broken minimum m 2 = (0.2 — 0.4)my. However, due to the cancellations, larger values 
are also possible. The remaining l/m 2 -divergence in eq. § is proportional to A, and 

1 /2 

its effect is damped even more by the fact that the l/m 2 -divergences come with the 
opposite sign. The dangerous term is -(l/167r 2 )(9^ 2 m^ /2 £ 3 / 4 )/4(2m 2 ) 1 /2 ; 

comparing 

it with the most important gauge-independent 2-loop terms, which are proportional 
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to g 2 rn^, we see that it is sufficient to have £ ^ (m^/m-p) 2 / 3 . The finite terms further 
suppress the divergent pieces at m 2 , 7^ 0, so that numerically their effect is confined to 
a very narrow region near m\ = 0, as is seen from fig. 2. 

Increasing £ still further, the loop expansion definitely breaks down. Namely, 
by just plotting the effective potential, one can see that the location of the broken 
minimum, (if), gets smaller for larger £. For very large £, this minimum gets so near 
the point m 2 = that it is in the region where the divergences are dominant. The 1- 
loop potential ceases to have any minimum at all in the region m 2 , > 0! To estimate the 
value of £ for which this happens, one can write all the ^-dependence in the equation 
dVi-i oop (ip, £ 3> \/g 2 )/dip = in terms of m\ to get (notice that m 2 < 0) 

e>W /2 \ 1/A Jrl'^uu = ^K-- 2 ) 1/2 +4A(3m 2 -2m 2 ) 1 / 2 -l^ m 2 . (9) 

By inspecting the functional form of the LHS and RHS of eq. |9| one sees that for large 
£, there is no solution with m 2 > (For smaller £, there are two solutions; the smaller 
of these is unphysical, see fig. 2). To get a numerical value, one can expand eq. [| to 
order m|, neglect even the m^-term for an upper bound on £, assume that A <C g 2 , 
and go somewhat below the temperature T . This leads to 



£ max ~ {g 3 /X 7/6 )(T/200V^) 1/3 , (10) 

with coupling constants scaled dimensionless by T. This means that for small Higgs 
masses and for temperatures near the point where the symmetric phase becomes 
metastable, £ can be quite large. 

Having estimated the range of validity of the 2-loop potential, we now study the 
convergence of the loop expansion. The object which hopefully converges well is not 
V((p,T,n) itself, but a renormalization group improved version of it. This means that 
one calculates dV((p, T, /j^/dtp, chooses \x for each ip so that large logarithms are killed, 
and then integrates to get an improved V(<p, T) f||, |14|| . In practice, one takes fi oc m T . 
One could kill the large logarithms directly in the potential V(ip, T, but in this case 
what used to be a /^-dependent vacuum part, becomes ^-dependent, and has to be 
taken into account |14| . 



Neither the effective potential at any temperature, nor its stationary points, are 
gauge-independent quantities || [L(| |llj] . However, the value of the effective potential 
at a stationary point is gauge-independent. This value is related to another gauge- 
independent quantity, ($^$). This is the quantity which we choose as our indicator of 
the convergence of perturbation theory: if we find at 2-loop level that ($^$) depends 
only weakly on the gauge parameter, we have reason to believe that the loop expansion 
converges rather well. 

For bare quantities, one gets the equation {&§)b = dVB((p m m{T), T)/dm 2 B . Un- 
fortunately, renormalization in a sense spoils this relation, since the (^-independent but 
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m 2 (/^-dependent counterterm in eq. ^ has to be removed. This also defines the renor- 
malization of (&$)b- As a result, (<3>t<3>) becomes /i-dependent, and thus physically 
ill-defined in the continuum theory. However, it can be related to a corresponding 
quantity on the lattice [ref. 0, eq. (13)], and thus given a physical meaning. Different 
values of \i just correspond to different lattices. To relate this //-dependence to the 
use of RG-improved perturbation theory, notice that the location of the minimum of 
V((p) is /i-independent (though ^-dependent) in the exact theory. Therefore, we use 
the RG-improved potential to find the best possible approximation for the location of 
the minimum, and then plug this value into the quantity dV(ip,T, fi)/dm 2 (fi) with a 
fixed /x, to get ($*$(//)). In practice, one cannot choose \i quite arbitrarily, since we 
are working at a finite order in the loop expansion. 



In fig. 3, the quantities (ip) and w2($t<I>) for m# = 80 GeV are shown at different 



temperatures below T c . One sees that y2($t$) is very accurately gauge-independent. 
We have also checked that there is less gauge-dependence at 2-loop level than at 1- 
loop level. One can think of two reasons for the remaining gauge-dependence. First, 
by comparing the most important terms in the effective potential, one sees that the 
effective loop-expansion parameter in the Landau gauge is g 2 T /mr{(p). Since at high 
temperatures (ip) gets smaller and T gets larger, the loop expansion converges worse, 
and there is more gauge dependence, as is observed in fig. 3. Second, the rather strong 
gauge dependence at large £ is caused by the fact that one is approaching the point 
where the loop expansion breaks down due to IR-divergences. 



In fig. 3, y2($t$) as calculated from the theory with Aq included, is also shown. 
The allowed range of £ is larger, as one can see from the same argument which led to 
eq. [K| There is a bit more gauge dependence in these curves, indicating that an error 
has been induced in constructing the Ag-integrated-out theory. 

Varying the Higgs mass, one can see that for smaller Higgs masses the allowed 
range of £ is much larger, as predicted by eq. [10]. The absolute value of the derivative 
d(&<&)/d£ is slightly larger for, say, m# = 35 GeV, than for m# = 80 GeV. 

Although it is not expected that the loop expansion converges very well in the 
symmetric phase, we have drawnQin fig. 4 the critical temperature T c and ((p c )/T c as a 
function of £. The gauge dependence is clearly stronger than for J2(&0>) deep in the 
broken phase. For large £ the loop expansion again breaks down. The reason is that at 
large £, T c is actually below the temperature T in eq. |3|. The logarithmic term makes 
m 2 {fi) positive for very small <p, since /2m is positive and \i is proportional to (p. When 
(p grows larger, m 2 (/i) becomes negative, but the term \ip 2 in m\ grows larger. The 
result is that for a large range of (p, m\ is very close to zero. Hence, the loop expansion 
is unreliable due to IR-divergences. 

To conclude, our results indicate that the loop expansion breaks down for too 
large £, but converges relatively well deep in the broken phase for smaller values. For 

2 At 1-loop level, a similar study was made in ref. JtI. See also ref. M. 
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more definite conclusions, one must compare perturbative results with lattice calcu- 
lations ||. Good convergence of the loop expansion, absence of IR-divergences, and 
sheer calculational simplicity, strongly support the use of Landau gauge. 
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